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The thermoelectric properties of intermetallic compounds with Ce or Yb ions are explained by the 
single-impurity Anderson model which takes into account the crystal-field splitting of the 4/ ground- 
state multiplet, and assumes a strong Coulomb repulsion which restricts the number of / electrons 
or / holes to n/ < 1 for Ce and n h f° le < 1 for Yb ions. Using the non-crossing approximation and 
imposing the charge neutrality constraint on the local scattering problem at each temperature and 
pressure, the excitation spectrum and the transport coefficients of the model are obtained. The 
thermopower calculated in such a way exhibits all the characteristic features observed in Ce and 
Yb intermetallics. Calculating the effect of pressure on various characteristic energy scales of the 
model, we obtain the (T,p) phase diagram which agrees with the experimental data on CeRu2Si2, 
CeCu2Si2, CePd2Si2, and similar compounds. The evolution of the thermopower and the electrical 
resistance as a function of temperature, pressure or doping is explained in terms of the crossovers 
between various fixed points of the model and the redistribution of the single-particle spectral weight 
within the Fermi window. 

PACS numbers: 75.30.Mb, 72.15.Jf, 62.50,+p, 75.30.Kz, 



I. INTRODUCTION 

The thermoelectric power, S(T), of intermetallic com- 
pounds with Cerium and Ytterbium ions exhibits some 
characteristic features which allow the classification of 
these compounds, into several distinct groups. 1,2,3 In the 
case of Cerium ions, the thermopower of the compounds 
belonging to the first group (type (a) systems) has a 
deep negative minimum at low temperaturesiiS^i 5 ^ and 
a high-temperature maximum, typically between 100 K 
and 300 K. At the maximum, S(T) could be cither posi- 
tive or negative, as shown in Fig. ^ At very low tempera- 
tures, the type (a) systems order magnetically or become 
superconducting. The compounds of the second group 
(type (b) systems) have a negative low-temperature min- 
imum and a positive high-temperature maximum but, in 
addition, the thermopower shows a smaller positive peak 
at lowest temperaturesi Vi 8 ^ 10 This second peak is some- 
times concealed by a low-temperature phase transition; 
for example, in CeCu2Si2 it becomes visible only in an ex- 
ternal magnetic field which suppresses the superconduct- 
ing transition^ and in CeRu2Ge2 it shows up when the 
external pressure suppresses the magnetic transition. 11 
The experimental evidence is now accumulating that 
the initial slope of the thermopower S(T)/T is positive 
for this class of (heavy fermion) materials, provided the 
measurements are performed at low enough temperature 
and with sufficient accuracyi§i2ii2 In the third group 
(type (c) systems), the low-temperature peak is well 
pronounced and shifted towards the high-temperature 
peak. The main difference with respect to the type 
(b) systems is that the sign-change of S(T) does not 
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Finally, in some cases (type (d) systems) 
the thermopower grows monotonically towards the high- 
temperature maximum, and the low-temperature struc- 
ture appears only as a shoulder on a broad peak, or is 



not resolved at alliVW 8 

The clue to these various types of behavior 
comes from the high-pressureSiiiii^iSSiSi and doping 
studies, 1 ! 14 ! 16 ! 24 ! 25 ! 26 ! 27 ! 28 which show that the ther- 
mopower of Cerium compounds changes continuously 
from type (a) to type (d). A typical example is provided 
by the S(T) of CeRu2Ge2, which is plotted in Fig^as a 
function of temperature, for various pressures^ At am- 
bient pressure, CeRu2Ge2 is a type (a) system with a 
magnetic ground state and negative thermopower below 
300 K. An increase of pressure leads to a thermopower 
with a small positive peak at low temperatures and an en- 
hanced peak at high temperatures. A further increase of 
pressure enhances both peaks, shifts the low-temperature 
peak towards the high-temperature one, and makes the 
thermopower at intermediate temperatures less negative. 
For large enough pressure, the sign-change does not oc- 
cur at all and for very high pressure the low-temperature 
peak merges with the high-temperature one, and trans- 
forms into a shoulder (see inset to Fig. The high- 
temperature peak grows continuously but its position re- 
mains more or less constant, as S(T) changes from type 
(a) to (c). Eventually, for pressures above 10 GPa, the 
S(T) assumes the (d) shape. Here, the initial slope of 
S(T) decreases and the position of the maximum shifts to 
higher temperatures, but its magnitude does not change 
as pressure increases. Similar behavior is also seen in 
the high-pressure data of, CeCu 2 Si 2 ri CeCu 2 Ge 2 | 20 i 22 or 
CePd 2 Si2i2 3 ' As regards doping, the substitutions which 
reduce the volume and make Ce ions less magnetic, trans- 
form S(T) from type (a) to type (b)j2£ from (b) to 
(c)r^ or from (a) to(c) l 25 i 26 i 27 while the substitutions 
which expand the volume and make the Ce more mag- 
netic, transform the thermopower from, say, type (d) to 
type (c) or from type (c) to type (b)j2£ This variation of 
shape is an indication that the local environment plays 
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FIG. 1: Temperature dependence of the thermoelectric power 
S(T) of CeRu2Ge2 for various pressures. Tk an d Ts label 
the centre of broad, pressure-induced maxima, related to the 
Kondo effect and the crystalline electric field, respectively. 
The inset shows S(T) data of CeRu2Ge2 in the non-magnetic 
phase. 



an important role in determining the magnetic character 
of Ce and Yb ions. Even at high temperatures, where 
each 4/ ion is an independent scatterer, the thermopower 
of a sample with a high concentration of 4/ ions cannot 
be obtained by rescaling the low-concentration data. 

The Ytterbium intermetallics can be classified using 
the mirror-image analogy with Cerium systems. This 
holds because the Yb ions fluctuate between 4/ 13 and 
4/ 14 , while the Ce ions fluctuate between 4/ 1 and if 
configurations, and the dynamics of a single / hole and 
a single / electron is the same. A well-defined local 
moment leads in Yb systems to the type (a) behav- 
ior, such that the thermopower has a negative minimum 
at high temperatures and a positive maximum at low 
temperatures)22i2S the size of the minimum is about the 
same as the size of the maximum. The thermopower of 
(b)-type Yb systems^ii 3 ^ 3 ^ mirrors the (b)-type Ce sys- 
tems. Here, one finds two negative minima separated 
by a small positive maximum. The type (c) Yb systems 
have a nonmonotonic thermopower with a large (nega- 
tive) minimum at high temperatures and a smaller one at 
low temperatures, but there is no sign-changei^Si^i* 3 ^ 3 . 
Finally, the thermopower with a single negative peak cen- 
tered around 100 Y^&SQiSLSL mirrors the type (d) behav- 
ior of Ce systems. The reduction of volume by pressure 



or dopin g 30 ! 33 stabilizes the magnetic 4/ configuration 
of Yb ions, and transforms S(T) from, say, type (b) to 
type (a), from (c) to (b), or from (c) to (a). 

The experimental results show that 4/ systems 
with similar thermopowers exhibit similarities in other 
thermodynamic^ and transport^LiSiiS properties, and 
there is an obvious correlation between the shape of S(T) 
and the magnetic character of the 4/ ions. The ther- 
mopower measurements provide a simple and sensitive 
tool for characterizing the magnetic state of a 4/ ion 
in a given metallic matrix: the shape of S(T) changes 
from the (a)-type in the case of magnetic Ce (Yb) ions 
with stable f 1 (4/ 13 ) configuration to the (d)-type for non- 
magnetic Ce (Yb) ions which fluctuate between the 4/ 1 
(4/ 3 ) and 4/> (4/ 14 ) configurations. 

We explain the thermoelectric properties of Ce and 
Yb ions in terms of a single-impurity Anderson model 
which takes into account the splitting of the 4/ states 
by the crystalline electric field (CF), and assumes an in- 
finitely large /-/ Coulomb repulsion, which restricts the 
number of / electrons or / holes to n/ < 1 for Ce and 
n hoie <- ^ f or yb. We assume that pressure changes the 
coupling and the relative occupation of the / and conduc- 
tion states, and impose the charge-neutrality constraint 
on the local scattering problem at each temperature and 
pressure. The total charge conservation provides a mini- 
mal self-consistency condition for a poor-man's treatment 
of pressure effects in stoichiometric compounds. The ex- 
citation spectrum of such a model in the vicinity of var- 
ious fixed points, the crossovers induced by temperature 
and pressure, and the corresponding effects on S(T), the 
number of /particles, n/(T), and the electrical resistance, 
p{T), are calculated by the non-crossing approximation 
(NCA) . The description of the stoichiometric compounds 
in terms of an impurity model is certainly inadequate at 
low temperatures where the / electrons become coherent. 
The errors due to such an approximation and the low- 
temperature errors inherent in the NCA calculations are 
discussed in detail at the end of Sec. [U] 

Our paper extends the long-standing theory of 
Coqblin^ii 4 ^ which described the high-temperature 
properties of Ce and Yb intermetallics by the Coqblin- 
Schrieffer (CS) model with CF splitting, and its more 
recent version 3 which improved the low-temperature cal- 
culations by rescaling the coupling constants. These pre- 
vious theories explained the main features of the temper- 
ature dependence of the thermopower^^ and the mag- 
netic susceptibility— but could not describe the pressure 
effects, because the CS model neglects charge fluctua- 
tions. Furthermore, the approximations used to solve the 
effective high- and low-temperature models cease to be 
valid at temperatures at which S(T) changes signj 3 . such 
that the shape of S(T) between the two maxima (min- 
ima) in Ce (Yb) systems could only be inferred from an 
interpolation. 

Here, we consider both the local spin and charge fluc- 
tuations, and provide the full description of the im- 
purity problem at various pressures and temperatures, 
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from above the CF temperature to below the Kondo 
temperature, including the intermediate regime where 
S(T) changes sign. Our results explain the shapes 
(a) to (d) of the thermopower, which are found in 
the systems like CeAl 3f 2 CeRu 2 Ge 2f ii^a CeCu 2 Si 2 ^i 
CeCu 2 Ge 2 , 20,22 or CePd 2 Si 2 , 23 and the 'mirror-image' 
shapes found in systems like YbNiSn, YnInAu 2 , YbSi 
or YbCu 2 Si 2 . 3 - We also explain the pressure data like 
on CeRu 2 Ge2ii or CeCu 2 Si 2 j2I and the chemical pres- 
sure data on Ce^Lai-^PdSn 27 and Ce a; Lai_ 2 ,Ru 2 Si 2 ^ or 
YbC^Sia^ 3 - 

The paper is organized as follows. In Sec. [H] we in- 
troduce the model, discuss its limitations, and describe 
the method of solution. In Sec. IIIII we provide the re- 
sults for the transport coefficients of Ce- and Yb-based 
intermetallics. In Sec. IIVI we discuss the effects of tem- 
perature and pressure on the spectral function, analyze 
the fixed-point behavior, and relate the shapes of the 
thermopower to the properties of elementary excitations. 
Sec. [V] gives the summary and the conclusions. 

II. THEORETICAL DESCRIPTION 

We model the intermetallic compounds by taking as 
many Ce or Yb ions per unit cell as required by the 
structure, but assuming that the scattering of conduc- 
tion electrons on a given 4/ ion does not depend on other 
4/ ions, except through the modification of the chemi- 
cal potential. In other words, we consider an effective 
impurity model which treats the 4/ states as scattering 
resonances rather than Bloch states but take into account 
the charge transfer due to the local scattering and adjust 
the chemical potential, fi, so as to maintain the overall 
charge neutrality of the compound. Such a description 
of the lattice problem applies at temperatures at which 
the mean free path of the conduction electrons is short 
and the scattering is incoherent. We consider mainly the 
Ce intermetallics and present only a few preliminary re- 
sults for Yb intermetallics. The Cerium ions are allowed 
to fluctuate between the 4/° and 4/ 1 configurations by 
exchanging electrons with the conduction band; the (av- 
erage) energy difference between the two configurations 
is \Ef\ and the hopping is characterized by the matrix 
element V. The Af configuration is excluded, i.e., an in- 
finitely strong Coulomb repulsion U between / electrons 
is assumed. The 4/ 1 configuration is represented by N 
crystal field levels: there are N — 1 excited states sepa- 
rated from the CF ground level by energies Aj -C \Ef\, 
where i — 1, . . . , N — 1. The local symmetry is taken 
into account by specifying the respective degeneracies of 
these levels, A/i. Thus, the low-energy excitations of Ce 
intermetallics are modeled by an effective single-impurity 
Anderson Hamiltoman^ 

H A = H b 

and ~t~ H imp ~l~ Umix-, (i) 

where H^and describes the conduction band, Hi„ lp de- 
scribes the CF states, and H m ix describes the transfer 



of electrons between 4/ and conduction states. In the 
absence of mixing, the conduction band is described by 
a semielliptical density of states, N(e), centered at 
and of half-width W, and the unrenormalized / states 
are represented by a set of delta functions at E*j and 
Ej — E® + Ai. The conduction states and the / states 
have a common chemical potential, which is taken as the 
origin of the energy axis. The properties of the model de- 
pend in an essential way on the CF splittings and on the 
coupling constant g = T/n\Ef\, where T = ttV 2 N(E®) 
measures the coupling strength between the / electrons 
and the conduction band, Ef = 53i=o -WiEf/Af, and 

Af = ^"q 1 ^, is the total degeneracy. We assume 
E° > 0, E f < 0, and T, A, < \E f \ ~ W, i.e., g <C 1. 
Since a single / hole is dynamically equivalent to a single 
/electron, we obtain the results for Yb ions by perform- 
ing the model calculations for a 8-fold degenerate / hole 
subject to the appropriate CF. 

The <? 1 limit of the Anderson model is controlled by 
several fixed points which are well understood. 43 In the 
case of a Ce ion with two CF levels split by A, the fixed 
point analysis can be summarized as follows. At small 
coupling, such that T < A <C \Ef\, we find nj(T) ~ 1 
and the model exhibits the Kondo effect. That is, all the 
physical properties depend only on the Kondo tempera- 
ture, To, which is uniquely determined by g, A, and the 
degeneracies of the CF states (for the NCA definition of 
To see Ref. — and Sec IIVI below). The low-temperature 
behavior is characterized by the Fermi liquid (FL) fixed 
point, which describes a singlet formed by an antifer- 
romagnetically coupled / electron and conduction elec- 
tron. An increase of temperature breaks the singlet and 
gives rise, at about To, to a transition to the local mo- 
ment (LM) fixed point, which describes a CF-split /state 
weakly coupled to the conduction band. The effective de- 
generacy of this / state is defined by the lowest CF level. 
For T > A, there is a further crossover to another LM 
fixed point, which describes the scattering of conduction 
electrons on a fully degenerate local moment. At higher 
coupling, such that T ~ 2 A <C \Ef\, the / charge is re- 
duced to 0.8 < nj(T) < 0.95, and the impurity still 
behaves as a Kondo ion, but the Kondo scale is much 
higher than in the case T < A < \Ej\. The two LM 
regimes are now close together and the crossover from 
the FL to the LM regime occurs at temperatures which 
are comparable to A. At very high coupling, such that 
2A < T < \Ef \ and n f {T) < 0.8, the /ions appear to be 
non-magnetic at all accessible temperatures due to the 
mixing of the 4j° and 4/ configurations. In this valence 
fluctuating (VF) regime, the behavior is non-universal 
and changes slightly, when the calculations are performed 
for different sets of parameters. Away from n/ ~ 1, more 
than one energy scale is needed to fully characterize the 
model. Other CF schemes pertinent to Ce and Yb ions 
in a different environment are characterized by similar 
fixed points. 

Our calculations show that the functional form of the 
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response functions changes at the crossover and that the 
g <C 1 limit of the Anderson model captures all the main 
features of the experimental results on Ce and Yb inter- 
metallics. To explain the pressure effects, which changes 
the thermopower of Ce and Yb systems in opposite ways, 
we assume that the exchange coupling g increases in Ce 
and diminishes in Yb compounds, as pressure increases. 
This difference arises because Ce fluctuates between 4/° 
and 4/ 1 , while Yb fluctuates between 4/ 14 and 4/ 13 con- 
figurations, so that the pressure-induced reduction of the 
number of electrons in the /-shell makes Ce ions less mag- 
netic and Yb ions more magnetic. 

In Ce intermetallics, there is a substantial overlap be- 
tween the / wave functions of Ce and those of the neigh- 
boring atoms, and we associate the pressure-induced 
increase of g ~ T/\Ef\ with an enhancement of the 
hybridization T. This enhances the Kondo tempera- 
ture, and pushes the system from the Kondo to the 
VF limit. In stoichiometric compounds, the pressure- 
induced reduction of nj is accompanied by the increase 
of n c , because the total charge of a given compound, 
ntot = n c + nf, is constant. The conservation of particles 
is enforced by adjusting fi, and since all the energies are 
measured with respect to /x, this amounts to shifting Ef 
and E c by some amount 5/i(T, T(p)). Thus, we describe 
the pressure effects for a given Ce compound by chang- 
ing F and keeping E c — Ef constant. The changes in the 
band-width and the CF splitting are neglected. 

In Yb intermetallics, the / states are more localized 
than in Ce systems, and we assume that the decrease in 
the radius of the 4/ shell as it loses charge at elevated 
pressure is sufficient to compensate for the increase in 
hybridization brought about by the reduction in unit cell 
volume. The reduction of g ~ T/\Ef \ in Yb compounds 
is achieved through an enhancement of the hole binding 
energy Ef as the neighbouring ions get closer to the rare 
earth, while T remains essentially constant. This reduces 
the Kondo temperature, and drives the system towards 
the Kondo limit. Since T is treated as a material-specific 
constant, we model the pressure effects in Yb systems by 
shifting Ef and solving for E c , so as to preserve n^°l e . 
This procedure shifts Ef and E c by different amounts and 
makes the separation E c — Ef pressure dependent. How- 
ever, when temperature is changed at constant pressure, 
the charge neutrality is enforced in the same way as for 
Ce compounds, by shifting the chemical potential with- 
out changing the separation E c ~ Ef. The changes in the 
bandwidth and the CF splitting are neglected. Describ- 
ing the pressure effects in such a way, we can calculate 
the response functions of the model for any value of the 
external parameters, and study the transitions between 
various fixed points. 

The electrical resistivity and the thermopower of the 
single-impurity Anderson model are obtained from the 
usual expressions^ 
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where L\\ and L\i are given by the static limits of 
the current-current and current-heat current correlation 
functions, respectively. In the absence of nonresonant 
scattering the vertex corrections vanish and the trans- 
port integrals can be written as^ 3 ^*^ 8 - 
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where Co is material-specific constant, f(uj) = 1/[1 + 
exp(w/fcs7 1 )] is the Fermi function, 1/t(lo) is the 
conduction-electron scattering ratej^ 3 . 
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A(lo) — T^Im Gf(u) ± i0 + ) is the /-electron spectral 
function, Gf{z) is the Green's function, and c is the con- 
centration of / ions. Eqs.QJ and JSJ show clearly that 
the sign and the magnitude of S{T) are determined by 
the spectral weight within the Fermi window (F window), 
i.e., by the shape of A{u>) for \lo\ < 2fcsT. The sign of 
S(T) is positive if the F window shows more states above 
than below the chemical potential, and is negative in the 
opposite case. The difficult part is to find Gf(uj ± i0 + ) 
and, here, we solve this problem by the NCA, following 
closely Refsi^ 3 ^, where all the technical details can be 
found. The main difference with respect to these NCA 
calculations is that we take c = 1 and enforce the overall 
charge neutrality. 

A detailed comparison with the experimental data 
shows that the transport coefficients obtained from the 
single-impurity Anderson model have all the hallmarks 
of the experiments, but discrepancies appear at low tem- 
peratures. This indicates the limitations of our approach 
which should be considered before presenting the NCA 
results. 

There are two main causes for the breakdown of the 
single-impurity model and the NCA calculations. First, 
at temperatures much below To the NCA spectral func- 
tion develops an unphysical spike, such that the resistiv- 
ity and the thermopower become artificially enhanced. 
This error becomes particularly severe at high pressure, 
because the characteristic scale To increases very rapidly 
with r and the non-analytic NCA spike appears at rather 
high temperatures. The unphysical enhancement of the 
low- frequency part of A(uj) reduces the integral for Ln, 
which is strongly underestimated at low-temperature. 
The integral L\i is less affected by this pathology, be- 
cause it has an additional w factor which removes the 
states within the Fermi window. Thus, the overall shape 
of S{T) ~ L? 2 CA /L? 1 CA seems to be qualitatively cor- 
rect, even though the low-temperature part of the curve 
has an unphysical enhancement. These difficulties are 
well known- 3 and relatively easy to resolve in the Kondo 
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limit, where the model has a unique low energy scale, To. 
We can find T in the LM regime, where the NCA is re- 
liable, and infer the low-temperature behavior from the 
universal power laws which hold in the FL regime. Com- 
bining the NCA results and FL theory, we can discuss 
the experimental data at all temperatures at which the 
single-ion approximation holds. 

A more serious problem is that, in stoichiometric com- 
pounds, the / electrons become coherent at low enough 
temperatures. This leads to a magnetic transition in 
CeRu 2 Ge 2 , the formation of a heavy FL in CeRu 2 Si 2 , su- 
perconductivity in CeCu 2 Si 2 , or some more complicated 
ground states. At high pressure, the coherence sets in 
at very high temperatures, as revealed by low values of 
the electrical resistance. The onset of coherence (like the 
NCA pathology) has its main impact on the low-energy 
states, giving l}^ Uce ^> jj^puny ^ gQ the impurity 
result badly overestimates the low-temperature electri- 
cal resistance. However, the considerations for the pe- 
riodic Anderson model^, or other models with on-site 
correlation^,, show that the integral for Ti 2 also con- 
tains an additional to factor which reduces the contri- 
bution of the low-energy coherent states to Li 2 , like in 
the single impurity case. Thus, our results for S(T) pro- 
vide a qualitative description of the experimental data at 
temperatures well below the onset of coherence, but the 
calculated values of S (T) around To are overestimated. 

We mention also that the analysis of the doping effects 
in terms of "chemical pressure" is not complete, because 
doping might give rise to a charge transfer or change 
the character of the ground state, and that the mirror- 
image analogy between Ce and Yb systems holds for the 
resonant scattering but may be lost in the presence of 
any additional scattering channel. Despite these draw- 
backs, the NCA solution of the Anderson model provides 
a surprisingly accurate description for a large body of ex- 
perimental data above the magnetic or superconducting 
transition temperature, and at low to moderate pressure. 



III. THERMOELECTRIC POWER RESULTS 

In this section we present the NCA results for the 
transport coefficients, describe the thermopower and the 
electrical resistance due to Ce ions in some detail, show 
the results for the / occupation, present a few prelim- 
inary results for the thermopower of Yb, and compare 
our results with the experimental data. 

Our strategy for Ce intermetallics is to illustrate the 
behavior of one particular compound as a function of 
temperature at various pressures. The compound is char- 
acterized at ambient pressure and high temperature by 
an initial parameter set {W, E®, N, T, E®, Aj, A/i}, where 
E® > and Ej<0 are measured with respect to [i 
and the high-temperature limit is defined by tempera- 
ture &bTa = Ajv-i- For given values of T and E° - Ej, 
we start the calculations at T = Ta and find dfi(T/^,T) 
such that the total charge is the same as the one obtained 
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FIG. 2: Thermopower of Ce ions calculated by the NCA for 
the CF splitting A = 0.07 eV is plotted as a function of 
temperature for several values of the hybridization strength 
r < 2A, as indicated in the figure. The bottom curve, T=0.06 
eV, corresponds to ambient pressure. 



for r = 0. (In the absence of coupling we have rif = 1 
and obtain n c by integrating the unperturbed density of 
states.) At high temperatures, the /-state is almost de- 
coupled from the conduction band, the renormalization of 
the parameters is small, and the numerics converge very 
fast. We then reduce the temperature, find the new shift 
Sfi(T, T) ensuring the charge conservation, and calculate 
the response functions for the resulting values of E c and 
Ef. This process is continued until the NCA equations 
break down at T <C Tq. To model the same system at dif- 
ferent pressure, we change T, find again <5/i(Ta,F), and 
repeat the same procedure as at ambient pressure for 
T < Ta- Note, n to t and E c — Ef are conserved at all 
temperatures and pressures. 

As a numerical example we consider a semielliptical 
conduction band of half-width W = 4 eV, centered at 
E® = 0.7 eV, and a / state split into a doublet and a 
quartet by the CF with A = 0.07 eV. We take n tot = 
n c + rt/=5.6301 electrons per ion (0.9383 electrons for 
each one of Af 'effective spin' channels). The transport 
coefficients are calculated for the hybridization strength 
changing from 0.06 eV to 0.20 eV, i.e., for T varying 
from r < A to T > 2A. The single-particle excitation 
spectra corresponding to these parameters are discussed 
in Sec. US 

Using the procedure outlined above we obtain for S(T) 
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FIG. 3: Thermopower calculated by the NCA for the CF split- 
ting A = 0.07 eV is plotted as a function of temperature for 
several large values of the hybridisation strength F, increasing 
from 0.13 eV to 0.2 eV. Fhe bottom curve, F=0.20 eV, is for 
the highest pressure. The inset shows the high-temperature 
maximum of S(T), plotted as a function of T. 



the results shown in Figs.[21and |3J The calculated curves 
exhibit all the shapes (a) to (d) found in the experimentas 
and give S(T) of the right magnitude, except at low tem- 
peratures where the calculated peak is too large with re- 
spect to the experimental data. As discussed already, the 
reason for this discrepancy is that the NCA overestimates 
the Fermi-level scattering rate for T <C Tq, and that we 
neglected the coherent scattering, which sets in at tem- 
peratures of the order of To. Thus, our low-temperature 
result for L\\ is artificially reduced, which makes S(T) 
too large. The sign and the topology of the S(T) curves 
do not seem to be affected by this error. 

For T < A, we have T < 5 K and n f > 0.95, and 
obtain S(T) with two well separated peaks, as shown 
by the r=60 and T=70 meV curves in Fig. The 
high-temperature peak is centered at Ts — Ta/2 and 
for our choice of parameters S max = S(Ts) < 0. The 
low-temperature maximum is at about To <C Tg and 
So = S(Tq) > 0. The thermopower between the two 
maxima is mainly negative. Since most of the type (a) 
and (b) systems order magnetically or become supercon- 
ducting above To, the low-temperature peak is not shown 
in Fig. |2| for T < 100 meV. A small increase of T (due to, 
say, an increase of pressure) reduces n / (see Fig. , en- 
hances Tq and S max , and expands the temperature range 
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FIG. 4: f-electron number, rif, calculated by the NCA for 
the CF splitting A = 0.07 eV is plotted as a function of 
temperature for several values of the hybridisation strength 
r. For the uppermost curve T=0.06 eV and then it increases 
in steps of 0.02 eV. At the bottom curve T=0.20 eV. 



in which S(T) is positive. Such a behavior, which is typ- 
ical of Kondo systems with small To, is in a qualitative 
agreement with the thermopower of the type (a) systems 
described in Sec[I] and with the data on CeRu2Ge2 n 
at low pressures (below 4 GPa) and above the ordering 
temperature, as shown in Fig^ 

For A < T < 2 A, we have T <150 K and nj > 0.8, 
and still obtain S(T) with the two maxima. But S max 
is now positive, the value of So is enhanced, the temper- 
ature interval in which S(T) < is reduced, and S(T) 
at the minimum is less negative than for smaller T. As 
we increase T (by increasing pressure), the sign-change 
of S(T) is removed, So and S m ax are further enhanced, 
but Ts is not changed. The two peaks are coming closer 
together and are merging eventually. These features are 
typical of Kondo systems with moderate To, say To > 10 
K, and are in a qualitative agreement with the data on 
the type (b) and (c) systems mentioned in Sec. U They 
are also shown by the CeRu2Ge2 data™ at intermediate 
pressures, (see the curves in Fig. ^ for pressure above 3.4 
GPa.) 

At T » 2A, the system enters the VF regime, S(T) 
acquires a single maximum at Ts, with a shoulder on the 
low-temperature side. Even this shoulder vanishes, when 
r is increased further, as shown in Fig. [3] For T > 2 A the 
thermopower is of the type (d), with a single peak which 
is much steeper on the high- than on the low-temperature 
side. This peak shifts to higher temperatures with in- 
creasing r (pressure) and S max saturates (see inset to 
Fig. EJ; the initial slope of S(T) decreases continuously 
with r. Such a behavior is in a qualitative agreement 
with the thermopower data on valence fluctuatorsi£iiLi£ 
and with the high-pressure data on CeRu2Ge2>ii' (cf. in- 
set to Fig. ^| and CeC^Siz^. However, a large dis- 
crepancy appears between theory and experiment at low 
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temperatures, because A(oj) is overestimated for u> ~ 0, 
which makes the NCA curves larger than the experimen- 
tal ones. A possible correction of the initial S(T)/T val- 
ues is discussed below. 

The /-electron number rif, calculated for the param- 
eters used in Figs. [21 and [31 is plotted in Fig. 21 as a 
function of temperature. The overall temperature depen- 
dence is rather slow, but two different types of behavior 
can still be seen. For T < 2 A <C —Ef, we find that rif is 
nearly independent of temperature and close to 1, which 
is typical of Kondo systems^ For T > 2A, we find that 
rif is less than 0.8 and nearly constant at low tempera- 
tures but at about T ~ A/3k B (271 K for A=0.07 eV) 
there is an increase followed by the saturation at high 
temperatures. Considered as a function of T (pressure), 
Tif(T) shows different behavior at high and low temper- 
atures. At high temperatures rif decreases uniformly as 
r increases. At low temperatures rif does not change 
much for r <C A and r » A, but drops rapidly around 
r ~ 2 A, indicating the crossover from the Kondo to the 
VF regime. 

The electrical resistivity, p m ag{T), obtained for small 
and intermediate values of T, is shown in Fig. The 
interesting feature is the high-temperature maximum, 
which appears for T < A, and correlates very well 
with the maximum in S(T). For temperatures below 
the maximum, p m ag(T) drops to a minimum and then 
rises logarithmically as Tq is approached. This mini- 
mum and the subsequent low-temperature upturn are of 
a purely electronic origin and appear in systems with 
small To and large CF splitting. In these systems, one 
can follow the evolution of the two peaks in p mag (T), 
and observe the disappearance of the minimum with the 
application of pressure. For example, in CeRu2Ge2, 11 
CePd 2 Ge 2f 2I CePd 2 Si 2! 2£ or CeCu 5 Au2a the minimum 
becomes more shallow, transforms into a shoulder, and 
vanishes at high enough pressure. As discussed already, 
the NCA overestimates the low-temperature scattering 
and distorts the relative magnitude of the high- and low- 
temperature peaks. In addition, the single-ion approx- 
imation always gives p m ag(T) which saturates at low 
temperature and cannot explain the low-temperature re- 
duction of p mag (T), which is seen in stoichiometric com- 
pounds below the onset of coherence. The electrical resis- 
tance of Ce-based Kondo systems at very high pressure, 
and the ambient pressure data of valence fluctuators, can- 
not be described by the NCA solution. In these systems, 
the scattering on / ions remains coherent up to rather 
high temperatures and the NCA solution is valid only 
above the high-temperature maximum and is not physi- 
cally relevant. 

To illustrate the situation in Yb intermctallics, we per- 
form a generic calculation for a scmielliptical (hole) con- 
duction band and four CF doublets. The NCA equations 
are solved for a single / hole (Af = 8). Following the 
procedure outlined in Sec. [HI we start the calculations 
at ambient pressure and temperature T<\, such that the / 




800 



FIG. 5: Electrical resistivity vs. temperature calculated by 
the NCA for the CF splitting A = 0.07 eV and for several 
values of the hybridization strength T, as indicated in the 



state is almost free, and calculate n^ e = n 1 l !: ole +n^ ole for 



figure. 



the initial parameter set {W, E°, N, T, Ej, A;, A/"i}. At 
lower temperatures we shift E c and Ef with respect to E® 
and E® by dp, so as to conserve n^°l e , and calculate the 
response functions for this new parameter set; this proce- 
dure does not change E c — Ef. At a higher pressure, we 
start again at T = Ta, change Ef so as to increase n^°* e 
and find E c which conserves n^ 1 ° 4 . Since T is not changed 
by this procedure, we now have E c — Ef ^ E® — Ef. For 
temperatures below Ta , the properties of the system are 
calculated by the same procedure as at ambient pressure, 
i.e., rit°l e is conserved by shifting E c and Et by the same 
amount 5 p. 

Taking W = 4 eV, E° c = 1.0 eV, T=0.08 eV, E° f = 
-0.8 eV, and three excited CF doublets at Ai=0.02 eV, 
A 2 =0.04 eV, A 3 =0.08 eV, respectively, we find n tot — 
6.444 at T = T A . The S(T) obtained for E f ranging 
from Ef = —0.6 eV to Ef = —1.0 eV, is shown in Fig. 
We recall, that an increase of pressure makes Ef more 
negative. For Ef = —0.6 eV, the thermopower shows a 
deep minimum, typical of Yb ions in the VF state. For 
Ef — —0.8 eV, S(T) develops a small maximum, which 
separates the high-temperature minimum at Ts and the 
low-temperature minimum at To. By making Ef more 
negative we shift To to lower values, as shown in Fig. EJby 
the curves obtained for Ef — —0.8 eV, Ef = —0.85 eV, 
and Ef = —0.9 eV, respectively. The low-temperature 
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FIG. 6: Thermopower due to Yb ions, obtained by the NCA 
for four CF doublets, is plotted as a function of temperature 
for various values of E®, i.e. for various pressures. Starting 
from the uppermost curve, we show S(T) for —E° = 1.00; 0.95; 
0.90; 0.85, 0.80, 0.70, and 0.60 eV, respectively. The bottom 
curve (Ef =-0.6 eV) corresponds to the lowest pressure. 



range in which S(T) is negative shrinks with pressure, in 
agreement with the experimental data. 30 However, for Et 
much below fx, the NCA calculations break down before 
this minimum is reached. As regards the value of S(T) at 
the maximum, it is negative at first but it becomes posi- 
tive as pressure increases, i.e., the thermopower changes 
from (c)-type to (b)-type. At very high pressures, such 
that the nj / 6 ~ 1, the thermopower is dominated at 
low temperatures by a large positive peak and at high 
temperatures by a negative minimum, which is typical of 
Yb-based systems with a small Kondo scale. The shape 
of S(T) is directly related to the magnetic character of 
Yb ions and our calculations explain the qualitative fea- 
tures of the thermopower of YbAu2, YbAu3r*i and YbSi, 
YnNi 2 Si 2 r& which are of the (a)-type, of Yb 2 Ir 3 Al 9! 3 ^ 
YbAuCu 4 ^i and YbNiSn^ which are of the (b)-type, of 
YbPdCu4^i and Yb2Rh 3 Al 9J 2£ which are of the (c)-type, 
and of YbAgCu4,YbPd 2 Cu 2 ri and YblnAu^ which are 
(d)-type. The pressure effects in YbSi^ and the chemi- 
cal pressure effects in YbCu 2 Si 2 33 are also in a qualitative 
agreement with our results. 

An estimate of the low-temperature properties of the 
single-impurity model can be obtained by combining the 
NCA results with the universal FL laws. The Sommerfeld 
expansion of transport coefficients in Eq. J3J gives 45 



urn — — --— = — — -sin — -rUf Z, o 

t^o k B T 3|e| Wo /; ' V ; 



where 7r 2 /cB/3|e|=283.5 fj,V/K and Z is the enhance- 
ment factor defined by the Fermi-level derivative of the 
/-electron self-energy, Z = [1 — dJ^/du] u=T _ . T/Z is 
related to the Kondo scale by a factor of order 1. Set- 
ting T = T/Z, and using for n f (T = Q) the NCA results 
of Fig. 01 we can estimate from Eq. @) the initial slope 
of the curves plotted in Figs. [21 and |21 The obtained 



values, which go from 32 /xV/K 2 for r=60 meV to 1 
^V/K 2 for r=130 meV and 0.1 fiV/K 2 for r=200 meV, 
are in the range reported recently for various Ce-based 
heavy fermions and valence fluctuators (see Table 1 in ref. 
— ). The initial slope of S(T) decreases as we move from 
the Kondo to the VF limit, in agreement with pressure 
experiments. 11 At higher temperatures, the non-linear 
corrections reduce S(T) and give rise to a maximum at 
Tq. These non- linear corrections are non- universal, and 
a large slope does not necessary translate into large S(T) 
at To. The corresponding calculations for the Yb curves 
plotted in Fig. show that an application of pressure 
enhances S(T)/T. 

The initial slope of S(T) can also be related to the 7 
coefficient of the specific heat, with the important result 
that the ratio S{T)/^T = (27r/|e|A/ )cot(7rn / /AA ) is in- 
dependent of Z. This expression and Eq.© are valid 
in the single-impurity regime, and it is not obvious that 
they would produce the correct results for stoichiometric 
compounds. However, the characteristic energy scales of 
a coherent FL (inferred from the experimental data on 
the initial slope of the thermopower, or the specific heat 
coefficient) do not seem to be much different from the 
single-ion scale To of the LM regime (inferred, say, from 
the peak in the thermopower or the Curie- Weiss temper- 
ature) . 



IV. DISCUSSION OF SPECTRAL PROPERTIES 

In this section we present the NCA results for A(u>), 
study the low-energy spectral features in the vicinity of 
various fixed points, discuss the changes induced by the 
crossovers, and explain the behavior of S(T) in terms 
of the redistribution of spectral weight within the Fermi 
window (F window). Only the Ce case is considered and 
it is assumed that pressure gives rise to an increase of the 
hybridization width T. The results obtained for T < A 
are shown in Fig. [7| where A(u>) is plotted as a function 
of frequency for several temperatures. At high tempera- 
tures, T ~ Ta, the spectral function has a broad charge- 
excitation peak somewhat above Ef and a narrower reso- 
nance of half- width A, centered below fi. This low-energy 
resonance is a many-body effect due to the hybridization 
of the conduction states with the 4/ states and is typi- 
cal of the exchange scattering on the full multiplet. In 
this temperature range, the F window shows more spec- 
tral weight below than above fi (see the middle panel in 
Fig. 01 and S{T) < 0. The magnetic susceptibility^ is 
Curie- Weiss like, with a very small Curie- Weiss temper- 
ature and a Curie constant which is close to the free Ce 
ion value. The maximum of S(T) at about Ts — T\/2 
is here negative, S max < 0, but a slight increase of T 
would make S max positive. All these features are typical 
of the LM fixed point corresponding to a fully degenerate 
/-state. At lower temperatures, T < Ts, the CF splits 
the many-body resonance into two peaks. The larger 
one grows below fi and the smaller one above fi (see the 
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FIG. 7: Spectral function A(ui), calculated for the hybridiza- 
tion strength T=0.06 eV and the CF splitting A=0.07 eV, 
plotted as a function of frequency for several temperatures. 
The solid, dashed, dashed-dotted and dotted curves corre- 
spond to T =2, 41, 209, and 670 K, respectively. The charge- 
excitation peak is visible in the upper panel. The middle 
panel shows the evolution of the CF and Kondo resonances 
with temperature. For T < A, the many-body resonance of 
half-width A is centered well below /j,. The F window has 
more states below than above /x and S(T) < 0. The lower 
panel shows the position of the Kondo resonance above /x. Its 
center defines To =1 K. For T < To the F window has more 
states above than below [i and S(T) > 0. 



middle panel in Fig. UJ. This asymmetry is enhanced 
as T is reduced, which is typical of the Anderson model 
with CF splittings; 4 — the increase of the low-energy spec- 
tral weight below \i gives rise to a large negative ther- 
mopower. A further reduction of temperature leads, for 
T <C 7a, to a rapid growth of an additional peak very 
close to fi, such that A(u>) acquires three pronounced low- 
energy peaks (see Fig. [7J. (The physical origin of these 
many-body resonances is explained in detail in Ref— 
The peak centered at wo "C A is the Kondo resonance 
and its appearance below T < 40 K marks the onset of 
the LM regime, which is due to scattering of conduction 
electrons on the lowest CF level. The two CF peaks cen- 
tered at about too ± A are outside the F window, and 
do not affect the low-temperature transport and ther- 
modynamics. Once the Kondo peak peak appears, the 
F window shows more spectral weight above than be- 
low /i and S(T) is positive, which is just the opposite to 
what one finds for T > TV The center of the Kondo 
resonance saturates at low temperatures at the energy 
luq > 0, which provides the NCA definition of the Kondo 
scale, ksTo — loq. In the symmetric Anderson model 
To is related to the width of the Kondo resonance but 
in the highly asymmetric case we are dealing with here, 
the current definition is more appropriate. The compari- 
son with numerical renormalization group (NRG) calcu- 
lations 50 shows that uio gives a reliable estimate of the 
Kondo temperature even for a doubly degenerate Ander- 
son model, and we assume that the NCA definition of To 
provides the correct Kondo scale of the CF-split single- 
impurity Anderson model as well. Because the Kondo 
resonance is asymmetric with respect to the u> = line 
and has more states above than below /x, the reduction 
of temperature enhances S(T) unltil it reaches, at To, 
the maximum value 5*o. A further temperature reduc- 
tion leaves the top of the Kondo resonance outside the F 
window, and the thermopower drops. However, most Ce 
and Yb system with very small To have a phase transi- 
tion above T , and to discuss the normal-state properties 
of (a)-type systems it is sufficient to consider the NCA 
solution for T > T . 

An increase of the coupling to T > A has a drastic ef- 
fect on A(oj), as illustrated in Fig. [SI where A(uS) is plot- 
ted for r = 0.12 eV. The charge-excitation peak is trans- 
formed into a broad background (see the upper panel in 
Fig. [BJ and the only prominent feature at T ~ Ta is 
the low-energy resonance of half- width A centered above 
fi. This low-energy resonance is due to the exchange 
scattering of conduction electrons on the full CF mul- 
tiplet, which gives rise to the maximum of S(T) in the 
LM regime. The F window (see lower panel in Fig. |SJ 
shows more spectral weight above than below /x, so that 
S(T) > 0. The reduction of temperature below T s re- 
moves some spectral weight above \x and brings addi- 
tional spectral weight below /x, which reduces S(T) and 
leads to a minimum— i A further reduction of tempera- 
ture leads to the rapid growth of the Kondo peak at too, 
and the CF peak at ujq + A, but the negative CF peak 
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FIG. 8: Spectral function A(lo), calculated for the hybridiza- 
tion strength F=0.12 eV and the CF splitting A=0.07 eV, 
plotted as a function of frequency for several temperatures. 
The solid, dashed, dashed-dotted and dotted curves corre- 
spond to T =2, 41, 209, and 670 K, respectively. The upper 
panel shows the overall features. The two many body reso- 
nances are resolved but the lower CF peak and the charge- 
excitation peak are absent. The lower panel shows the evolu- 
tion of low-energy resonances with temperature. For T > A, 
there is more spectral weight above than below at all tem- 
peratures and S(T) is always positive. 



does not develop. That is, an increase of pressure re- 
moves the lower CF peak, and shifts the Kondo and the 
upper CF peak to higher energies, without changing their 
separation A. The F window shows more spectral weight 
above than below fi, so that S(T) is positive and grows 
as temperature is lowered. The maximum So is reached 
at To when the Kondo resonance is fully developed. The 
characteristic energy scale is defined again by the posi- 
tion of the Kondo peak, fc^To = ujq, which can now be 
quite large. For T <Tq the F window becomes narrower 
than the Kondo resonance and S(T) drops below Sq. For 
T < T , where the FL behavior is expected^ the NCA 
gives A(u>) with an unphysical spike at [i, which makes 
Pmag{T) and S(T) much larger than the exact result. 
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FIG. 9: Spectral function A(u>), calculated for the hybridiza- 
tion strength T=0.140 eV and the CF splitting A = 0.07 eV, 
plotted as a function of frequency for several temperatures. 
The solid, dashed, dashed-dotted and dotted curves corre- 
spond to T=2, 41, 209, and 670 K, respectively. For F > A, 
the Kondo resonance is reduced to a small hump above the / 
level. 



However, once To is obtained from the NCA calculations, 
the low-temperature transport can be inferred from the 
universal power laws which hold in the FL regime, as 
discussed in the previous section. 

A further increase of T shifts the Kondo and the CF 
peaks to higher energies, and changes their relative spec- 
tral weight, as shown in Fig. where the low- frequency 
part of A(lo) is shown for T = 2A. The Kondo scale 
is still defined by the center of the Kondo peak, even 
though it is now reduced to a hump on the low-energy 
side of a large peak centered at uq + A. The unphysical 
NCA spike at lu — can be seen at lowest temperatures. 
The thermopower is positive at all temperatures and has 
only a shoulder below Ts. A quantitative comparison 
between To, defined by the position of the Kondo reso- 
nance, and the position of the Kondo anomaly in S(T) 
becomes difficult. 

Finally, for T > 2 A, we find A(cu) with a single broad 
peak centred at Ef > 0, as shown in Fig. ^| The CF 
excitations are now absent, which is typical for the An- 
derson model in the vicinity of the VF fixed point. The 
relevant energy scale at low temperature is defined as 
UbTq = Ef, and shows an almost linear dependence on 
F. The thermopower is always positive and grows mono- 
tonically from small values at low temperatures towards a 
high-temperature maximum at Ts. The unphysical spike 
at lu = appears at higher temperatures and is more 
pronounced than for small T, as illustrated in Fig. ^| 
The initial slope of S(T) obtained from the NCA result 
for transport coefficients is very much overestimated with 
respect to the FL result based on Eq.©. An increase of 
temperature above Ts modifies the excitation spectrum 
on an energy scale of the order of Ef and reduces S(T). 
The values of Ef and Ts do not seem to be related in any 
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FIG. 10: Spectral function A(u>), calculated for the hybridiza- 
tion strength T=0.20 eV and the CF splitting A = 0.07 eV, 
plotted as a function of frequency for several temperatures. 
The solid, dashed, dashed-dotted and dotted curves corre- 
spond to T=2, 41, 209, and 670 K, respectively. For T > A, 
the Kondo resonance is absent. 
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FIG. 11: The temperature Ts, defined by the high- 
temperature maximum of the thermopower, and the char- 
acteristic scale To, defined by the peaks in the spectral func- 
tion, are plotted as functions of the hybridization strength V. 
The three regions, where To changes its functional form, are 
clearly seen. The inset reproduces the experimental data for 
CeRu 2 Ge 2 and CeCu 2 Ge 2 p±i 



simple way, and the temperature-induced breakdown of 
the local FL is not indicated by any FL scale. An increase 
of r (pressure) changes A(oj) and enhances Ts but does 
not lead to further enhancement of S max - The pressure 
dependence of Ef and Ts looks different. In other words, 
the description of the VF regime requires more than one 
energy scale. 

The characteristic scales of the Anderson model ob- 
tained for various values of T provide the (T,p) phase 
diagram of the system plotted in Fig. ^2 The To line is 
defined by the position of the Kondo peak and the Tg- 
line by S max - At small T (low pressure) Fig. ITT1 indicates 
two crossovers. The one around To is between the FL 
and the LM regime defined by the lowest CF level. The 
one around Ts is the crossover to the high-temperature 
LM regime, as defined by a full CF multiplet. When 
these LM regimes are well separated, the crossover be- 
tween them is accompanied by a minimum and, possi- 
bly, a sign-change of S(T). At intermediate pressures, 
the two LM regimes are too close for the sign-change to 
occur, and the crossover is indicated only by a shallow 
minimum or just a shoulder of S(T). Here, the rela- 
tionship between the low-temperature maximum of S(T) 
and the center of the Kondo peak can only be given as 
an order-of-magnitude estimate. At very high pressures, 
the system is in the VF regime and the crossover from 
a universal low-temperature FL phase to a non-universal 
high-temperature phase takes place at T ~ Tg. This 
crossover is not defined by the FL scale Ef, which is 
very large, but by a much smaller scale Ts- 



V. CONCLUSIONS AND SUMMARY 

We have applied the single-impurity Anderson model 
to the investigation of the temperature and pressure 
dependence of the thermopower of Ce and Yb inter- 
metallics, and found that the crossovers between various 
fixed points explain the seemingly complicated tempera- 
ture dependence of S(T). 

The basic assumption of our approach is that for a 
given concentration of rare-earth ions, and above some 
coherence temperature, the system is in an "effective 
impurity" regime. In that "impurity limit", we treat 
the rare-earth ions as independent scattering centers and 
solve the ensuing single-impurity model by the NCA. Our 
calculations impose the charge neutrality constraint on 
the local scattering problem, which provides a minimal 
self-consistency condition for describing the pressure ef- 
fects. The excitation spectrum obtained in such a way 
is very sensitive to the changes in the coupling constant 
g = T/n\Ef\, and an increase of the hybridization or a 
shift of the / state, which modifies the magnetic state 
of the / ion, has a huge effect on the spectral function. 
Since the excitation spectrum is related to the transport 
coefficients by Kubo formulas, the thermopower changes 
rapidly as a function of temperature, pressure or doping. 
The shapes (a) to (d) assumed by S(T) in various regions 
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of the parameter space follow straightforwardly from the 
redistribution of the single-particle spectral weight within 
the Fermi window. 

The NCA solution of the single-impurity Anderson 
model breaks down for T < Tq. However, the excitations 
above the ground state are determined by the general FL 
laws, which can be used to find the initial slope of S(T). 
Estimating the Kondo scale and the number of / elec- 
trons by the NCA calculations, we find that an increase 
of pressure reduces (enhances) the low-temperature value 
of S{T)/T in Ce (Yb) compounds. Combining the FL 
and the NCA results, we can obtain S(T) in the full 
temperature range at any pressure. 

Our results explain the thermopower of most inter- 
metallic Ce compounds mentioned in Section [i] The 
multi-peaked S(T), which characterizes the type (a) 
and (b) systems, is obtained for Kondo ions with small 
To and well-defined CF resonances. A type (c) ther- 
mopower, which has weakly resolved peaks or just a 
broad hump, is obtained for Kondo ions with large To and 
partially overlapping CF resonances. A single-peaked 
S(T) of the (d)-type is obtained for the VF ions which 
do not show any CF splitting of the single-particle ex- 
citations. If we assume that pressure increases T, and 
reduces nj, our results account for the changes of To, 
So, S(T)/T, and S max , observed in pressure experi- 
ments on Ce compoundsi 2 ! 11 ! 19 ! 20 ! 21 ! 22 ! 23 The strong dop- 
ing dependance 1 of S(T) is explained as a chemical pres- 
sure effect, which changes the / occupation. In our local 
model the " effective concentration" of / electrons is de- 
termined self-consistently, and the thermopower of con- 
centrated and dilute systems cannot be related by simple 
scaling. (Just like the high pressure data cannot be de- 
scribed in terms of rescaled ambient pressure data.) The 
relevant energy scales obtained for different values of T 
agree with the (T,p) phase diagram of CeRu2Ge 2 and 
CeCu2Ge2, inferred from the high- pressure transport and 
thermodynamic data^i 

For Yb compounds, we argue that pressure or chemi- 
cal pressure mainly affect the position of the crystal field 
levels relative to the band center, while leaving the hy- 
bridization width unchanged. By shifting Ef, so as to 



increase the number of holes, we obtain the thermopower 
of the type (a) to (d), in agreement with the experimental 
data. The qualitative features seen in Yb intermetallics 
at various pressures or chemical pressures are captured 
as -w-elljSLSii^iiSiiS but different compounds require dif- 
ferent initial parameters, and a detailed analysis is yet to 
be done. 

In summary, the normal-state properties of stoichio- 
metric compounds with Ce and Yb ions seem to be very 
well described by the local model which takes into ac- 
count spin and charge fluctuations. The classification 
of the thermopower data follows straightforwardly from 
the fixed point analysis of the single-impurity Anderson 
model with the CF splitting. A rich variety of shapes 
assumed by S(T) at various pressures or doping dramat- 
ically illustrates the effect of the local environment on 
the response of the single rare earth ion. The nature 
of the ground state, and the fact that the ground state 
of many Ce and Yb systems can be changed with pres- 
sure or doping, do not seem to affect the thermopower 
in the normal state. We take this as an indication that 
the normal state properties of Ce and Yb compounds are 
dominated by the local dynamics. 

The main problem with our poor man's treatment of 
ordered compounds is that it neglects the coherent scat- 
tering which sets in at low enough temperatures, and 
reduces the thermopower below the values predicted by 
the local Fermi liquid theory. We can argue that these 
effects do not change the qualitative features of the ther- 
mopower but the proper answer cannot be obtained with- 
out solving the lattice model. 
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